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Chiral magnetic effect in the absence of Weyl node 
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The nodal points in a Weyl semimetal are generally considered as the causes of the chiral anomaly 
and the chiral magnetic effect (CME). Employing a linear-response analysis of a two-band lattice 
model, we show that the Weyl nodes and thus the chirality are not required for the CME, while 
they remain crucial for the chiral anomaly. Similar to the anomalous Hall effect, the CME results 
directly from the Berry curvature of energy bands, even when there is no monopole source from 
the Weyl nodes. Therefore, the phenomenon of the CME could be observed in a wider class of 
materials. Motivated by this result, we suggest that the nodeless CME may appear in three- 
dimensional quantum anomalous Hall insulators, but after they become metallic due to the band 
deformation caused by inversion symmetry breaking. 
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I. INTRODUCTION 

Materials with topologically-nontrivial electronic 
structure have recently been under intensive investiga¬ 
tion. A particularly interesting state of matter is the 
three-dimensional Weyl semimetal,whose band struc¬ 
ture contains isolated band-touching points, called Weyl 
nodes. Such nodal points behave as sources/drains of the 
Berry flux and carry nonzero monopole charges Q, which 
is defined by the integral of the Berry curvature Ok over 
the surface enclosing the node, Q = (l/27r)^dSk ■ Ok. 
The Weyl nodes are protected topologically against 
perturbations. They would disappear only if two nodes 
with opposite monopole charges merge and annihilate 
with each other. Recently, TaAs and NbAs have been 
experimentally confirmed to be Weyl semimetals.®® 
This progress paves the way for exploring novel effects 
in these materials. 

Due to the nontrivial topology in momentum space, 
such nodal materials exhibit a wide variety of unusual 
electromagnetic responses.^ In a Weyl semimetal, elec¬ 
trons near each Weyl node can be assigned a chirality by 
the monopole charge of that node. Applying a pair of 
non-orthogonal electric and magnetic helds, the charges 
can be transported between two Weyl nodes with oppo¬ 
site chiralities. Therefore, the number of electrons with 
a definite chirality is no longer conserved, showing the 
so-called chiral anomaly.^ By using a semiclassical anal¬ 
ysis, one can show that the anomalous source term in 
the continuity equation of chiral charges is proportional 
to the monopole charge.l^ Thus this exotic phenomenon 
becomes more manifest for nodes with larger Q’s. This 
anomaly is predicted to give an enhanced negative mag¬ 
netoresistance when the applie d elec tric and magnetic 
fields are parallel to each other.l^^l^ Such a prediction 
has just been confir med by experiments in the Weyl 
semimetal TaAs.li®^^ 

Besides chiral anomaly, such nodal materials may show 
chiral magnetic effect (CME ) whe n the energies of pairs 
of Weyl nodes are different 7H 2ni This effect gives a 



FIG. 1: (color online). Various types of semimetals, (a) Two 
conical bands that tonch at a nodal point. The Fermi surface 
is a point, (b) Tilted cones (or type-II Weyl semimetaP^. 
The Fermi surfaces enclose an electron pocket and a hole 
pocket, (c) Overlapped bands that do not touch with each 
other. Dashed lines denote the location of the chemical po¬ 
tential. 


dissipationless electric current J flowing along an applied 
magnetic field B, J = — oB. Employing a low-energy 
effective theory with unbounded linear dispersion, the 
CME coefficient a is shown to be proportional to the 
energy separation between a pair of Weyl nodes. Since 
the CME can be related with the chiral an omaly through 
the energy balance of chirality generation ,1121211 one may 
expect that, when there is no Weyl node such that chiral 
fermions become ill-defined, neither the chiral anomaly 
nor the CME could exist. 

However, most of the early investigations are based 
on effective models with linear dispersions around Weyl 
nodes. In a previous work, we find that, when going be¬ 
yond the linear regime such that the concept of chirality 
may no longer be appropriate, the CME can still exist .1^ 
Our analysis shows that this effect is better understood 
in terms of the Berry curvature, rather than the chirality. 

In this paper, we show that, in sharp contrast to the 
chiral anomaly, Weyl nodes are in fact not vital for the 
existence of the CME. As shown in Fig. there are 
several ways to realize a semimetal. In Figs. Sa) and 
(b), when the nodes carry nonzero monopole charges, 
chiral anomaly and related transport phenomena are ex¬ 
pected to emerge. In contrast, since there is no nodal 
point in Fig. [^c), chirality becomes meaningless and the 
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monopole charge density (i.e., the divergence of Berry 
curvature) vanishes throughout the whole Brillouin zone. 
Therefore, there is no chiral anomaly in this nodeless 
case. However, we find that for clean samples the CME 
can exist even in this case with no Weyl node. 

This conclusion may not be too surprising because, for 
both cases in Figs, [^b) and (c), energy bands are par¬ 
tially filled with finite Fermi surfaces, and the system is 
metallic in nature. Besides, in the semiclassical analy- 
sis,l^ there exists anomalous velocity induced by a Berry 
curvature: = —k x fik, where /ik is the force expe¬ 

rienced by the electron. Under an external B field, the 
Lorentz force gives an anomalous velocity (and thus cur¬ 
rent): Va ~ I (v • J7k)B, where v is the group velocity of 
the Bloch electron. While this provides only a heuristic 
picture [the complete expression is presented in Eq. ([^], 
it does show that the Berry curvature could drive a cur¬ 
rent along the direction of the B field, no matter whether 
the Berry curvature emanates from a monopole or not. 

Our observation is illustrated by using a two-band lat¬ 
tice model. With suitable parameters, cases (a), (b), and 
(c) in Fig. can all be realized. We find that, when 
two nodes merge in momentum space [see Fig. [^a)] and 
the resulting node has zero total monopole charge, the 
CME would appear as long as there are electron and 
hole pockets at the Fermi level. That is, the energy sepa¬ 
ration between two nodes and nonzero monopole charge 
of the nodes are not necessary for the CME. This case is 
reminiscent of the type-II Weyl semimetal in Fig. 0 b)- 
Furthermore, for some parameters, the point degener¬ 
acy between energy bands can be lifted, so that there 
is no Weyl node. However, the CME still would appear 
when the energy bands are partially filled [see Fig.j^b)]. 
Our results challenge the wisdom based on the studies of 
linearized models, in which the chiral magnetic current 
arises from the presence of Weyl nodes and is propor¬ 
tional to the energy separation between Weyl nodes. 

This paper is organized as follows: In Sec. H, we in¬ 
troduce the two-band model under consideration. In 
Sec. HI, the dependence of the anomalous Hall conduc¬ 
tivity and the CME coefficient on various parameters is 
studied numerically. The result of this work is summa¬ 
rized in Sec. IV. 


II. TWO-BAND LATTICE MODEL FOR 
DOUBLE-WEYL SEMIMETAL 

Here we consider the case of double-Weyl semimet¬ 
als for illustration. Besides the usual linear Weyl nodes 
carrying monopole charges Q = ±1, there are nodal 
points with nonlinear dispersions and higher monopole 
charges.l^^l^ Such nodes can be protected by crystallo¬ 
graphic point group symmetries. One example of the so- 
called double-Weyl semimetal is HgCr 2 Se 4 , which con¬ 
tains nodes with Q — ±2. The spectrum around each 
node disperses quadratically in two directions. This ma¬ 
terial has recently been confirmed by a transport ex¬ 


periment to be a semimetal.^ It is predicted that, in 
a double-Weyl semimetal, the quantum anomalous Hall 
(QAH) conductivity, the coefficient of the chiral anomaly, 
and the number of Fermi arcs for surface states could all 
be doubled. The effect of electron interaction near such 
a quadratic node could also be very different from that 
of a linear node.^^^H^ 

To study the CME of double-Weyl semimetals, we start 
with a two-band lattice Hamiltonian motivated by the 
compound HgCr 2 Se 4 ,l^ 

i7o(k) = (cos kx — cos ky) -|- sin kx sin kyU^ 

+ {m — cos kx — cos ky — cos kz) cr^ , ( 1 ) 

where (t“ {a = x, y, z) are the three Pauli matrices 
and k is the Bloch wavevector. This model breaks time- 
reversal symmetry, but has a combined C'4-rotation and 
Mj-mirror symmetry. The bulk energy gap of Eq. Q 
closes and a pair of double-Weyl nodes emerges if 1 < 
|m| < 3. When 1 < to < 3, the nodes locate at kg = 
(0, 0, ± cos“^(to — 2)); while for —3 < to < —1, they 
locate at kj = (tt, tt, 7r± cos“^(to-I- 2 )) . Notice that the 
C4 symmetry protects the double-Weyl points, while the 
Mz symmetry requires that the two double-Weyl points 
have equal energy. 

Since we are interested in studying the CME, a Mz- 
breaking term, Usinfez, is added to split the energy of 
the two double-Weyl nodes. Besides, in order to split 
a double-Weyl node into two single-Weyl nodes, a C 4 - 
breaking term, acr^, is added. Therefore, the lattice 
model under consideration becomes 

H (k) = HoCk) + a < 7 ^ + ti sin kz 

= d(k) ■ (T + ti sin kz , ( 2 ) 

where we have written H in the standard d-vector nota¬ 
tion. The energy spectrum are invariant under the com¬ 
bined transformation of a —> —a, and kx ky. There¬ 
fore, it is sufficient to study the case with a > 0 . 

The phase diagram in the parameter space of (to, a) is 
shown in Fig. The ti term only tilts the band struc¬ 
ture, but does not change the number of nodes in a phase. 
When a = 0, as mentioned above, the system belongs to 
the semimetallic phase with a pair of double-Weyl nodes 
for 1 < |to| <3. If 1 < TO < 3, and a 7 ^ 0, then each of 
the double-Weyl node would split to two nodes along the 
fc^-direction. The resulting phase has four linear Weyl 
nodes within the colored parallelogram in Fig. Simi¬ 
larly, if—3 < TO < — 1 , then a positive a would split a 
double-Weyl node to two nodes along the fcj^-direction, 
and we have another Weyl semimetal phase with four 
nodes. These two four-node phases would overlap to form 
a phase with eight nodes at larger a’s. 

On the other hand, if |to| < 1 and a is small (the 
upright triangular region in Fig. [^, then there exists no 
node. When ti = 0 and the chemical potential ^ = 0, 
the system becomes a three-dimensional QAH insulator 
with a quantized Hall conductivity jgffj = 2e^//i (here 
the lattice constant is set to be unity) .^3321 The coefficient 
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FIG. 2: (color online). Phase diagram of the Hamiltonian in 
Eq. § for a > 0. Here, a circled number (in blue) shows 
the number of Weyl nodes in a particular phase. Doubly cir¬ 
cled numbers indicate that they are double-Weyl nodes (see 
Appendix A for more details). This occurs when a = 0 and 
1 < \m\ < 3 (two red line segments). Outside the colored re¬ 
gions, the system is a trivial phase without node. The vertical 
dashed line indicates the path calculated in the right panels 
of Fig. |4] 


of two can be traced back to the fact that the monopole 
charge of a double-Weyl node is two times larger. 

By choosing suitable parameters, both cases of tilted 
cones and overlapped bands in Figs, [^b) and (c) can be 
realized in the present model. For example, when a = 0 
and m = 1, two (would-be) double-Weyl nodes merge at 
kg = (0,0,7r). Notice that the merged node carries no 
monopole charge. For a nonzero ti with a fixed chemical 
potential /i = 0, electron and hole pockets appear at 
the Fermi energy, as shown in Fig. |^a). This resembles 
the type-II Weyl semimetal discussed in Ref. [23l On the 
other hand, there is no node when the parameters (m, a) 
lie within the upright triangular region in Fig.[^ Turning 
on a nonzero ti, the energy bands are distorted so that 
the chemical potential fi = 0 intersects with both bands. 
A typical example is displayed in Fig. ib). 


III. ANOMALOUS HALL EFFECT AND 
CHIRAL MAGNETIC EFFECT 


Within the linear-response theory, both the Hall con¬ 
ductivity and the CME coefficient can be calculated from 
the retarded current-current correlation function at finite 
frequency and wavevector, followed by an appropriate 
limiting process.!^ 

For a generic two-band model, the Hall conductivity 
a’’^ in the relation is 

given as 

where i,j,£ are in the cyclic order of x,y,z. Here /t(k) 
is the Fermi-Dirac distribution function for band t[= ±) 
at some temperature T. The Berry curvatures are 
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FIG. 3: (color online). Band energies as Inactions of kz {k^ = 
ky). The parameters are (a) m = 1 (a single merged node 
with zero topological charge) and (b) m = 0.8 (non-touching 
energy bands). Here a = 0 and ti = 1. 


given by the formula,^ 
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k,± 
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4d3(k) 


d(k) 


dd(k) dd(k) 
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dki 


dki 


■ ( 4 ) 


Since the Weyl nodes in our model are connected by Dirac 
strings along the z-axis, the only non-vanishing compo¬ 
nents are = —<jh- The rest of the compo¬ 

nents are zero, which has been confirmed by numerical 
calculations. 

For a two-band model, following the steps described in 
Ref. [22 one can show that the CME coefficient for the 
chiral magnetic current JqMp = — under an applied 
B held along direction i ip4l 


e 
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d^k 

(27r)3 
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Vk,+ + Vk,- 
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—t (i(k) 


^^k,t /t(k) 

t—zk ^ 

Vk,t • dft 

2 dEt 


( 5 ) 


Here Vk,t = (l/h)VkEt(k) are the group velocities, and 
Et{k) are the band energies. Because the integrands of a® 
along different directions differ only by a term containing 
the values of a* do not differ much. They are 
often of the same order and show similar dependence on 
parameters. (This is true also for the model considered 
in Ref. |22l) For simplicity, only its averaged value, a = 
{ax -\- tty + Uz)/i-, is presented below. 

The numerical results of ur and a for the two-band 
model in Eq. ([^ are displayed in Fig.|^ We start with the 
0 = 0 case. When E — 0, the system is nodeless and be¬ 
longs to the QAH insulating phase with a quantized Hall 
conductivity an = —2e^lh if |m| < 1.1^ As m increases 
from 1 to 3, the two double-Weyl nodes approach each 
other, and thus the magnitude of <jh decreases to zero, 
when they annihilate each other at m = 3. When ti ^ 0, 
such that the valence band is not completely filled, or 
is no longer quantized even when |to| < 1. These results 
are summarized in Fig. |^a).l^ 

On the other hand, as seen from Fig. |^b), the CME 
coefficient a vanishes as long as ti = 0. This is expected. 
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FIG. 4: (color online). Hall conductivity an (in units of fh) 
and CME coefficient a (in units of jh^). Left panels: a = 0 
and m £ [0,3] (along the positive a-axis in Fig. [^. Right 
panels: m = 0.5 and a £ [0,2] (along the vertical dashed 
line in Fig. [^. The parameter ti varies from 0 to 1 with 
step 0.25. Temperature T = 10“^ and the chemical potential 
fj, — 0. Phase boundaries at H = 0 are denoted by the vertical 
dotted lines. The calculations are done with 600^ lattice sites. 


because the system is either in the insulating phase 
(|m| < 1) or the nodes has no energy shift (1 < m < 3). 
Notably, when ti 0, a can be nonzero even if the sys¬ 
tem lies in the nodeless region of |m| < 1. The finite 
value of a comes from the effect of the Berry curvatures 
of the electron and hole pockets [see Fig.j^b)]. The value 
of a drops to zero at a smaller m (a larger energy gap) as 
long as the chemical potential no longer intersects with 
the energy bands. Our results shows that the Weyl node 
is not required for the CME. 

Now we consider the effect of the C' 4 -breaking term 
with a yf 0. In Figs, ^(c) and (d), one follows the dashed 
line of TO = 0.5 in Fig. that starts in the phase with 
non-touching bands when a = 0 , enters the four-node 
phase when a = 0.5, and reaches the eight-node phase 
when a = 1.5. At ti = 0, there are kinks in the slope 
of the Oh curve at phase boundaries in Fig. I^c) .1221 The 
curves become smooth when ti yf 0. Interestingly, as seen 
in Fig.gd), a behaves smoothly and non-monotonically 
when ti yf 0. Due to the same reason mentioned in the 
previous paragraph, a = 0 as long as ti =0. On the two 
ends of a = 0 and a = 2 , an and a can be non-zero at 
non-zero ti, because the valence band is not completely 
filled. Most importantly, again the CME coefficient is 
obviously not zero in the nodeless region when a < 0.5. 

Lately, SrSi 2 is found to be a double-Weyl semimetal 
with pairs of nodes located at different energies.l^Sl Thus 
it could be a natural candidate for testing the CME. 


IV. CONCLUSION AND DISCUSSION 


ters, it exists even in the absence of Weyl node and the 
chirality becomes ill-dehned. Our analysis shows that the 
CME comes directly from the distribution of the Berry 
curvature in energy bands. Whether there are monopoles 
(sources/drains of the Berry flux) in the Brillouin zone is 
not crucial. 

It worths emphasizing that the existence of the CME 
in a type-II Weyl semimetal or a nodeless phase is not 
restricted to the systems containing double-Weyl nodes. 
Same phenomena are found to exist (not shown here) 
also in a lattice model (an extended version of the one in 
Ref. I 22 II with only linear Weyl nodes. Because the CME 
coefficient is odd under space inversion, only the mate¬ 
rials without inversion symmetry can exhibit the CME. 
Inspired by the present work, a candidate for realizing 
the nodeless CME could be a three-dimensional QAH in¬ 
sulator, after its energy bands being distorted by, e.g., 
inversion symmetry breaking, so that the valence (con¬ 
duction) band is not completely filled (empty). 

Recently, the chiral anomaly in Weyl semi-metal 
has been confirmed indirectly through the measure¬ 
ment of the negative longitudinal magnetoresistence 
(NLMR). b^ft^ lln the experimental setup, one applies elec¬ 
tric and magnetic fields to a Weyl semi-metal. Because 
of the chiral anomaly, charges are pumped between Weyl 
nodes with opposite chiralities at a rate proportional to 
E • B, which leads to a difference in chemical potentials 
between these Weyl nodes. Since the CME coefficient 
is proportional to such difference of chemical potentials, 
the electric conductivity is thus proportional to re¬ 
sulting in the NLMR. However, in this paper, we are 
studying the CME in an equilibrium state without an 
electric field. That is, the electric current is driven only 
by a magnetic held. Thus there is no connection between 
the CME here in an equilibrium state and the NLMR. In 
the case of nodeless CME, even if an additional electric 
held is applied, since there is no node, no chiral anomaly 
and associated charge pumping, so one does not expect 
to see the NLMR either. 

Finally, we emphasize that the present linear-response 
analysis is done with a two-band model in a clean and 
inhnite system. Real samples always have disorders and 
are hnite in size. An important issue is how would these 
factors change the value of the CME coefficient. For ex¬ 
ample, the spin Hall conductivity for a clean Rashba sys¬ 
tem is ftrst predicted to be non-zero, but later found to 
vanish with the inclusion of disorders.^ Similar qualita¬ 
tive alteration to the CME due to realistic factors cannot 
be completely ruled out. In any case, the conclusion in 
this paper should still hold in the conservative range of 
LOT ^ 1 , where lo is the driving frequency and r is the 
relaxation time. Generalization to more realistic systems 
remains to be explored in future investigations. 


Even though the CME is originated from the theory of 
chiral fermions in high-energy physics, in condensed mat- 
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Here one sums over the k-points with d(k(^)) = do, which 
is fixed, assuming that the Jacobian \ddi/dkj\ 7 ^ 0 at 
these k-points. 

For example, for a = 0 and m = 2 in our model, there 
are nodes at kg = (0, 0, ± 7 r/ 2 ). Expanding the momen¬ 
tum around the node, k = kj -f q, one has 


Appendix A: Calculating the monopole charge of a 
Weyl node 


d(q) 



2 J 


(A2) 


The energy dispersion near a Weyl node determines its 
monopole charge. As discussed in Ref. [371 the monopole 
charge Q (or the winding number of the mapping / ; k —>■ 
d) of a given Weyl node can be calculated by a simple 
formula: 


N 

Q^^sgn 

e=i 


ddi 


dki 


k^kCC) 


(Al) 


and the Jacobian is 

^ =T{ql + ql)- (A3) 

Choosing dp = (1/2, 0,0), then there are two q points, 
q(i) = (0,1, 0) and q^^^ = (0, —1,0). They contribute to 
Q = =f 2 in total. 
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